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'l«^', Defect electrons (holes) play an important role in most technologically important complex oxides. 

Qs ' In this contribution we present the first detailed characterization of localized hole states in such 

O^ \ materials. Our investigations employ advanced embedded-cluster calculations which consistently 

include electron correlations and defect-induced lattice relaxations. This is necessary in order to 

account for the variety of possible hole-state manifestations. 
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I. INTRODUCTION 



o 

C/3 ] The basic structural building units of solid oxides are MOg metal-oxygen octahedra and, possibly, additional MO4 

tetrahedra. The complete crystal structure is built up of corner-, face- and/or edge-sharing connections of these 
structural elements. Additional cations (A) can be incorporated at interstitial lattice sites. These are increasingly 
G formed, the more open-structured the whole network of MOn units appears to be. Generally, open crystal structures 

• _ imply high formal M-cation charge states (referring to formal 0^~ anions) leading to mixed ionic-covalent (or semi- 
d . ionic) material properties. The M-cations are in most cases transition-metal (TM) ions. Examples of such complex 
H [ materials are given by AMO3 perovskite-structured oxides such as barium titanate (BaTiOa). Perovskite oxides 
are frequently ferroelectric and possess important electrooptic applications based on the photorefractive effect. Also 
high-Tc oxides resemble the perovskite structure. 
Q . Charge carriers, either valence-band (VB) holes or conduction-band (CB) electrons, are created by doping with 

CJ ■ impurities, annealing treatments or by light-induced charge-transfer excitations which, for example, take place dur- 
ing photorefractive processes in the appropriate oxides. Combined optical-absorption- and electron-spin-resonance 
,__! ■ measurements (Photo-ESR) Q proved that the created and afterwards trapped holes influence the light-induced 
^ t charge-transfer reactions in photorefractive BaTiOa. These holes, which are either para- or diamagnetic, probably 
Cn ' induce the observed sublinear dependence of photoconductivity on the light intensity (e.g. |g]). A further example 
v^ highlighting the relevance of holes refers to high-Tc oxides: Pairs of doped holes give rise to superconductivity in 

^»^ these materials. The possible pairing scenarios are still a matter of active debate. 

^^ . The present work initiates a systematic and detailed characterization of localized hole states in complex oxides. 

f — ' Our investigations employ real-space embedded-cluster calculations, which consistently combine electron correlations 

On . and defect-induced lattice relaxations. This procedure is indispensable in order to predict the richness of possible hole 

states. Here, we demonstrate simulations for trapped holes in BaTiOs, but many results, ranging from stabilization 

of cationic charge states to the formation of bipolarons, can be extrapolated to other oxides. 

S " The trapping of holes at acceptor defects leads to a localization of hole states, which is further aided by defect- 

1-^ , induced lattice distortions. In near-insulating high-T^ oxides the antiferromagnetic order supports hole localizations. 

r^ ' Principally, there is 'on-acceptor' and 'off-acceptor' hole trapping. In the first case the trapped hole localizes at 

O . the acceptor site (thereby increasing the acceptor charge state by one positive unit, i.e. M+" -I- h* — > M+("+^)), 

J-^ ■ and in the second case the hole remains at the oxygen ligands. Further, the off-acceptor case allows to classify hole 

^ , states according to their degree of delocalization, i.e. complete localization at exactly one ligand anion, intermediate 

Lj localization at two neighboring oxygen ions (Vk centers), and delocalization over more than two oxygen ligands. 

''^ Which hole-localization type is favored depends on the ionicity of the host system, but also on the incorporation 

d ' site and electronic structure of the acceptor. For example, at Ba-site acceptors localized off-acceptor holes are likely 

to be preferred over delocalized species due to large interionic separations [0. Further, the simultaneous binding 

of two holes at one acceptor facilitates the formation of small intersite bipolarons representing negative-C/ centers. 

Such species could explain the diamagnetic hole centers in photorefractive oxides |l||; intersite bipolarons have also 

been considered as possible superconductivity-carrying bosons in high-Tc oxides (see [y, for details). Earlier work 

on bipolarons in oxide materials dealt exclusively with electron type species. We quote in this context the pioneering 

work of P. W. Anderson ||5|] on the formation of negative-C/ centers in amorphous solids. Schlenker [|6[ investigated 

small electron type bipolarons in titanium oxides. 



II. METHODS 

Our investigations employ real-space embedded-cluster calculations (ECC), in which the crystal is devided into 
a central and quantum chemically treated defect cluster and into the embedding lattice which is modeled more 
approximately. The ECC discussed in this publication deal with the trapping of holes at Ti-site acceptors in BaTiOa. 
Our calculations are based on 21-atom clusters MOeBagTig centered around an acceptor-oxygen octahedron MOe, 
see also 0. All simulations assume the perfect-crystal structure of cubic BaTiOs. This simplification may be justified 
due to the observation that all possible ferroelectric distortions of the material are small compared with usual defect- 
induced lattice relaxations. 

The molecular-orbital ansatz for the central MOg complex employs Gaussian type basis functions with split- valence 
quality for the acceptor metal impurity and its oxygen ligands; the oxygen basis set is further augmented by polarizing 
d-functions. Bare effective core potentials [B| are used to model the localizing ion-size effects of the outer Ba- and 
Ti-cations. The ab initio level of our calculations covers Hartree-Fock (HE) theory, M0ller-Plesset perturbation theory 
to second order (MP2) and density functional theory (DFT). All quantum cluster simulations have been performed 
on the basis of the quantum-chemical CADPAC code M. Within DET two choices are used to approximate the 
exchange-correlation functional, i.e. the local spin density ansatz of Vosko, Wilk and Nusair m^ (VWN-LSDA) 
and the generalized-gradient-approximation functional involving the Becke exchange term |11[] and the correlation 
functional derived by Lee, Yang and Parr ||l2|. This functional, abbreviated as 'BLYP', improves on LDA. 

The modeling of the embedding lattice and of cluster-lattice interactions employs a shell-model pair-potential 
description (see also Refs. [[7| |l3|jl4[| ) . The shell-model parameters have been taken from the earlier work of Lewis and 
Catlow [|l5| . It is noted that the shell model takes electronic polarization contributions of the embedding lattice ions 
reasonably well into account. This feature is particularly important with respect to the highly polarizable oxygen 
anions. Eigure |l| displays the embedding scheme. As in any classical Mott-Littleton defect calculation the (embedding) 
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FIG. 1. Visualization of embedded-cluster calculations. The quantum cluster is embedded in a classically treated lattice 
consisting of the regions I, Ila and lib. The arrows indicate the various existing interactions between the crystal regions. 



crystal lattice is divided into three regions: region I which is explicitly equilibrated according to the underlying pair 
potentials, the interface region Ila and region lib. Region II(a and b) is treated as a polarizable continuum using the 
harmonic approximation for the region-II self-energy. Whereas the interactions between the cluster and region I, on 
the one hand, and region Ila, on the other hand, are included explicitly, all region-lib species feel only the effective 
defect charge of the central defect cluster. Details of this well-known description are given in Ref. |l6|. All lattice 
calculations employ the CASCADE code 0. 

The total energy of the complete model crystal is given by the expression: 



E(p,Rc,RiJ 



p,Hc,Hr 



EEnv(I^r) -f E '' {p, Re, Rr)} - 



(1) 



p, Re and Rr denote the cluster electron density, the coordinates of the cluster nuclei and the positions of the shell- 
model species in the embedding environment, respectively. Egiust represents the quantum mechanically calculated 
cluster self-energy, E|^^ the Mott-Littleton type shell-model energy of the embedding environment and E^"* is the 
cluster-lattice interaction energy. The required minimization in Eq.(Q) is performed using the programs CADPAC 
(variation of p with fixed Rr) and CASCADE (variation of Rr with fixed Re and p). In order to do cluster geometry 
optimizations which are consistent with a predetermined embedding crystal lattice (represented by a point-charge 
field) an additional program has been written which updates the total cluster energies and gradients, as calculated by 
any quantum chemical program such as CADPAC, by adding the appropriate short-range pair-potential contributions 
due to the interactions between cluster ions and embedding lattice species. With these updates the program carries 
through the cluster geometry optimization using a variable-metric (quasi-Newton) minimization algorithm. The 
total minimization procedure is completed, when the main cycle consisting of the alternating lattice- and cluster- 
equlibration subcycles converges. It is finally noted, that we have neglected the full cluster multipole consistency 
during the embedding-lattice relaxation step (see |lq] , for details), since, in our simulations, we did not observe any 
appreciable relaxation and energy changes due to cluster multipole contributions deviating from formal-charge models. 
The present approach proves to be the only successful one which can consistently integrate the local electronic 
structure and large-scale lattice-distortion effects. In particular, we find the combination of DFT and of our embedding 
scheme to be sufficiently accurate and flexible in order to realistically account for the variety of possible hole states. 
We stress that cutting off the electronic structure beyond the cluster boundary does not principally invalidate our 
results. This is due to the fact that embedded-cluster calculations formally refer to a localized-orbital point of view 
(e.g. see |19 3] for related discussions) and not to the extended one-electron eigenstates of the crystal. Moreover, all 



hole states considered in this contribution may be assumed to be localized within the cluster region due to the presence 
of the acceptor defect and due to pronounced lattice deformations (small-polaron effect). Further indications towards 
a reliability of our present results are provided by our earlier calulations of optical charge-transfer (CT) transitions 
at transition-metal cations 0, which compare favorably with experimental data. The CT states may be considered 
as excited states of trapped holes. 

III. RESULTS AND DISCUSSION 

In the following discussion, in each case we start with the ionic HF model and add on correlations afterwards to 
demonstrate important implications. 

First, we discuss single holes trapped at TM cations. Within HF theory, on-acceptor holes are favorable for 
divalent TM cations including Cu^+, which possess filled electron levels above the VB edge. At trivalent Ti-site 
TM acceptors, on the other hand, localized (i.e. symmetry-breaking) off-acceptor holes are highly preferred (Table 
||). Here, we observe the hole localization in 7r-type oxygen 2p-orbitals (perpendicular to the oxygen-acceptor axis). 

TABLE I. Total HF and MP2 embedded-cluster energies for different hole states including lattice relaxations. Energies are 
renormalized to localized off-acceptor states. Oq{-k, a): localized off-acceptor hole in oxygen 2p n- or a-type orbitals. (0(xy))~: 
off-acceptor hole delocalized over the four oxygen ligands in the xy-plane. 2S -I- 1 denotes the total spin multiplicity of the 
cluster. Note that Fe^+ and Cr^+ favor high spin (2Sfo -1-1 = 6, 2Scr -1-1 = 4), and Rh''+ prefers low spin (2SRh 4-1 = 2). 

Defect 2S + 1 AE^„Hi^(eV) AEZM^V) 

Fe^+ - Oo(7r) 7 Ol) OO 

Fe^+-0o(7r) 5 0.0 



Ti 

Fei+ 5 +2.4 +0.1 



Fe:^+-Oo(cr) 5 +1.0 

Fe^l - (0(xy))- 7 +2.4 +0.5 



'^Ti 



Rh:^7-Oo(7r) 2 0.0 0.0 

Rh^+ - (0(xy))- 2 +2.7 +0.1 

Rhi| 2 +0.4 -2.4 






Cr:^t-0o(7r) 5 0.0 0.0 

Cr*+ 2 +1.15 -2.6 

This symmetry-broken state is stabilized by defect-induced lattice distortions and, in this way, it can acquire physical 
significance. The relaxation-induced energy gain is typically 1 eV, which emphasizes the importance of such terms. 



The occurrence of symmetry-broken solutions may be considered as electronic polaron self-trapping. Generally, hole 
localization in cr-type orbitals (parallel to the oxygen-acceptor axis) is less favorable, but deviations may occur due 
to the actual electronic acceptor structure. 

Interestingly, both delocalized off-acceptor holes and the formation of tetravalent TM ions are unfavorable within 
HF. For example, delocalized off-acceptor holes are by 2 — 3 eV discriminated against their localized counterparts. 
Similar results have been obtained in previous HF simulations on trapped single holes in MgO P^,pi[ and NiO [p2| . 
Also these calculations predict localized off-acceptor holes to be most favorable. In fact such holes are likely the most 
favorable off-acceptor solutions in highly ionic oxides like MgO, but this behavior may prevail in any HF calculation, 
thus simulating artificially ionic situations. 

In BaTiOs the effects of electron correlation are most significant for trivalent TM acceptors due to the increased 
covalent charge transfer from oxygen onto the TM ions. This has important consequences: The on-acceptor hole 
localization is enhanced, and, simultaneously, possible off-acceptor holes increasingly delocalize. In contrast, at most 
divalent TM cations including Cu^+ on-acceptor holes remain stable. This behavior differs from La2Cu04 where, 
confirmed by our preliminary ECC, doped holes localize at the oxygen sites. The different effective copper charges in 
both oxides can explain such deviations. Further work on La2Cu04 is in progress and will be published in the near 
future. 

Our MP2 calculations (Table ||) illustrate that electron correlation tends to favor on-acceptor holes at trivalent 
TM acceptors. Rh^+ and Cr"'+ become preferred over M — O^ centers, but Fe''+ remains 0.1 eV less favorable than 
Fe^"*" — 0~. Within DFT even Fe^+ is stable, and for all trivalent TM acceptors the most favorable off-acceptor holes 
are delocalized at this level. 

Considering off-acceptor holes, correlation mostly affects the delocalized hole states, but localized states to a much 
lesser extent. These off- acceptor- hole results demonstrate the pronounced interplay between orbital relaxations and 
electron correlation, i.e. the localized hole states with significant HF orbital relaxations receive less correlation- 
energy gain than delocalized hole states. However, different from HF theory the energy separations between localized 
and delocalized off-acceptor holes remain within some tenths eV. This result indicates that a modest increase of 
ionicity could favor localized over delocalized off-acceptor holes. The lanthanide contraction of 5d TM cations might 
accomplish such changes leading to observable localized off-acceptor holes, if these are sufficiently stable against the 
formation of on-acceptor holes. ESR data suggest that localized off- acceptor holes bound to ESR-silent Pt-acceptors 
exist in BaTiOs [0. The charge state of platinum could not be assessed experimentally, but Pt^+ with a completely 
filled 5d{t2g) subshell is a possible candidate. 

In summary, electron correlation essentially stabilizes the high-valent charge states of various TM cations in BaTiOs. 



This mechanism resembles the Haldane- Anderson scenario |23 in semiconductors. In highly ionic materials, in 
contrast, the charge-state stabilization has been ascribed to lattice relaxations p3. The present discussion confirms 
that BaTiOa is semi-ionic. 

In comparison to most TM acceptors the hole localization properties can differ at ionic acceptor cations due to the 
lack of filled acceptor levels above the VB edge and of covalent charge transfer. Table H exemplifies the situation 
for Al'^^. Obviously, on-acceptor hole states leading to Al**"^ are highly unfavorable. ESR investigations furnished 

TABLE II. Hole formations at Alrpt. Total embedded-cluster energies are renormalized to the localized off-acceptor 7r-type 
hole state. The results include lattice relaxation. 

Defect AEli^ljeV) AE?^a(eV) AEg.^rP(eV) 

aItI -05(7r) ao oo oo 

A4+ - (deloc.(7r))- +2.5 4-0.4 -0.8 

Al^t - (deloc.(cr))- +3.7 +1.9 

A4+ - (Vk(7r)) +1.0 0.0 -0.4 

Al^t - (Vk(cr)) +2.8 +1.6 

evidence in favor of Vk hole centers (i.e. one hole delocalized over two neighboring ligands) ||25| , but there are no such 
indications for other acceptors in BaTiOa. 

As expected, HF favors localized off-acceptor 7r-type holes. Next favored is the Vk 7r-type hole distribution (Fig||), 
but it is by 1 eV less favorable. The situation is even worse for completely delocalized states. Again electron 
correlation supports hole delocalization. The correlated calculations reported in Table || confirm that all hole states 
are of comparable energy. Most favorable are the delocalized state and the formation of Vk centers. However, there is 
no clear indication in favor of the latter species. But the results suggest a delicate balance between the local electronic 
structure and defect-induced lattice distortions: Increasing the acceptor ionicity would lead to localized off-acceptor 





FIG. 2. Hole spin density plot of n- and a-type Vk centers trapped at Alrpj . The occupation of oxygen 2p orbitals can be 
clearly seen. The acceptor is at (x,y) = (0,0) in both cases. All xy-coordinates in the figures are given in atomic units. The 
interatomic separation between the oxygen partners forming the Vk defect is 2.3 A and 2.6 A for localization in tv- and cr-type 
oxygen orbitals, respectively. The perfect-lattice separation is 2.8 A. 



holes, whereas a reduced ionicity implies hole delocalization. In an intermediate small 'window' the formation of Vk 
centers might be most favorable. This also explains why Vk centers are rarely observed in oxides. Experimentally the 
Al-cation fits into this window, but also the present embedded-cluster model is close to this situation. Interestingly, 
Vk centers have been also proposed in AI2O3 [ p6[ . 

Due to reduced bonding cr-type Vk centers are by 1.6-1.8 eV less favorable than 7r-type ones. All ion displacements 
are smaller in the cr-type modification. We stress that these results refer to Al'^+ acceptors, but it is recalled that 
holes doped into the Cu-0 planes of high-Tc oxides are assumed to populate cr-type oxygen orbitals |g^, which, 
following our preliminary calculations, is related to the occupation of copper 3dx'^-y2 states and to the presence 
of strong electron correlations. Our present results receive particular importance, since Vk centers (upon addition 
of a second hole) bear resemblance to small intersite hole bipolarons. Thus, we may expect the same relaxation 
behavior for hypothetical 7r-type and a-type hole bipolarons in high-Tc oxides. Due to moderate lattice distortions 
the latter species would be more appropriate to coherent motions. Speculatively, cr-type bipolarons could facilitate 
superconductivity, thereby supporting the scenario of Alexandrov and Mott 0]. Certainly, our speculations deserve 
further detailed investigations. 

We finally review trapped hole bipolarons in BaTiOs; the details of the calculations are discussed extensively in 
Refs. 13; 14]. Ti-site Mg^+ is considered as the actual acceptor cation. Two geometrical hole-acceptor configurations 
are compared: the bipolaron with two holes on neighboring oxygen ions (i.e. literally a Vk center with one additional 
hole) and a linear complex O^ — MgrpT — O^. The presence of two holes introduces a natural driving force towards 
localization, of which the linear complex corresponds to minimizing the inter-hole Coulomb repulsion. Both hole 
complexes have been studied under different conditions: (1) HF treatment of the cluster employing a rigid and 
perfectly structured crystal lattice. Only the actual O^ partners are allowed to relax. (2) HF description of the 
cluster including lattice relaxation. (3) Correlated calculations (MP2 and DFT). To save computer capacities further 
geometry optimizations have been performed at the DFT level (LSDA and BLYP) only. 

We first survey the perfect-lattice HF simulations. Electronic interactions between the O^ ions and nearby crystal 
ions favor the spin triplet over the singlet state. Due to its antisymmetrical charge distribution the triplet-state 
bipolaron benefits most from this interaction. Increasing the bipolaron bond length rapidly disturbs the bipolaron 
type states due to hole-state delocalizations. But energetically most favorable is the linear complex. This expected 
result reflects the disfavor of delocalized holes within HF. The linear configuration with either spin is 1 eV more 
favorable than the triplet-state bipolaron. 

Next we consider the effect of lattice deformations. Most importantly, our ECC demonstrate that defect-induced 
lattice relaxation and electronic correlations stabilize hole bipolarons in BaTiOs. Lattice relaxations increase the 
localization of bipolaron states and favor the spin-singlet state. Thus, such relaxations enable the possibility of 
embedded Oj" molecules analogous to isoelectronic F2 dimers. Noticeably, both species are unstable within HF (e.g. 
see 10). The ultimate stability is established by correlation-induced bonding terms. Binding energies of bipolarons 
can be estimated as the energy difference E(BP) — E(0^ — Mg — 0~), since, due to shell-model estimates, the 



competitive linear complex is only weakly bound in BaTiOs. We obtain -0.41 eV (MP2), -1.13 eV (DFT-BLYP) 
and -2.1 eV (DFT-LSDA). Calculated bond lengths are: 1.48 A (LSDA) and 1.55 A (BLYP). As is frequently observed 
p8[ , LSDA overestimates bonding. The accurate binding energy is enclosed by the BLYP- and MP2 values. 

We believe that paired holes are of general importance in any oxide. Possible differences will refer to the occupation 
of TT- or CT-type oxygen 2p orbitals giving different bonding strengths. Finally, we emphasize that the singlet-triplet 
splitting of bipolaron states, i.e. the 'spin gap', depends on the strength of lattice relaxations. Whereas in perfect 
lattices the triplet state is preferred, the singlet state becomes most favorable upon lattice distortions. For 7r-type 
bipolarons in BaTiOs, due to our MP2 calculations, the spin gap becomes 1.50 eV in this case. We note that 
this value is essentially intrinsic to the hole bipolarons in BaTiOs, and does not depend on the specific nature of 
the magnesium acceptor ion. In particular, the bipolaron states are not contaminated by magnesium orbitals; we, 
therefore, expect a similar spin gap for possibly existing isolated bipolarons in BaTiOs. In their bipolaron theory 
of high-T(7 superconductivity Mott and Alexandrov suggested a spin gap of a few tens meV Q . Such a small value 
would indicate that, compared with BaTiOs, lattice relaxations are present but less developed in high-Tc oxides. 
This observation is important if coherent motion of bipolarons is to be required. 
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